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ABSTRACT 

We present a multiwavelength study of the Galactic H Ii region Sh2-297, located in Canis Major 
OBI complex. Optical spectroscopic observations are used to constrain the spectral type of ionizing 
star HD 53623 as BOV. The classical nature of this H ii region is affirmed by the low values of electron 
density and emission measure, which are calculated to be 756 cm~^ and 9.15x10^ cm~^ pc using 
the radio continuum observations at 610 and 1280 MHz, and VLA archival data at 1420 MHz. To 
understand local star formation, we identified the young stellar object (YSO) candidates in a region 

of area ~ 7.5 x7.5 centered on Sh2-297 using grism slitless spectroscopy (to identify the Ha emission 
line stars), and near infrared (NIR) observations. NIR YSO candidates are further classified into 
various evolutionary stages using color-color (CC) and color-magnitude (CM) diagrams, giving 50 red 
sources {H — K > 0.6) and 26 Class Il-likc sources. The mass and age range of the YSOs are estimated 
to be ~ 0.1-2 Mq and 0.5-2 Myr using optical (V/V-I) and NIR (J/J-H) CM diagrams. The mean 
age of the YSOs is found to be ~ 1 Myr, which is of the order of dynamical age of 1.07 Myr of the 
H II region. Using the estimated range of visual extinction (1.1—25 mag) from literature and NIR 
data for the region, spectral energy distribution (SED) models have been implemented for selected 
YSOs which show masses and ages to be consistent with estimated values. The spatial distribution 
of YSOs shows an evolutionary sequence, suggesting triggered star formation in the region. The star 
formation seems to have propagated from the ionizing star towards the cold dark cloud LDN1657A 
located west of Sh2-297. 

Subject headings: dust, extinction - H ii regions - ISM: individual objects(Sh2-297) - Infrared: ISM - 
Radio continuum: ISM - Stars: formation 



1. INTRODUCTION 

The study of massive stars (> 8 Mq) is one of the 
most prolific areas of research in present day astronomy, 
a fact which can be attributed to the not-so-well under- 



stood processes and a range of theories (see Zinnecker & 

I |kshitiz@tifr.res.in| I 



Yorke|2007 and references therein) . Massive star forma- 
tion, though a rare event in itself, can have a significant 
impact on its natal environment through processes such 
as strong stellar wind and supernovae explosions, which 
transfer large energy and momentum to the surrounding 
environment. Features such as bright-rimmed clouds and 
H II regions (owing to their large output of UV photons) 
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are hallmar ks of massive sta rs in a region. 
( [Sharpli 



Sh2-297 
Figure [T]) 

star HD 



■less 
classical 



1959) is an optically visible (see 
Galactic H ii region (a20oo 
12° 19 00 ), ionized by a massive 



, (52000 

53623, and bounded by a cold, dark cloud 
LDN1657A to the west. It is associated with the re- 
flection nebula IRAS 07029-1215 and is a part of Canis 
Major (CMa) OBI association. In the literature, the dis^ 
tance estimates for S h2-297 range from 1 to 1.4 kpc (Felli| 



k Harten|[l98T 
In the p resent s 



Bica et al.]|2003 [Forbrich et al.| |2C 
udy, we "Have adopted the value ot 1.1 



kpc from Bica et al. (2003). Ruprecht (19661 estabhshed 
the boundaries of CMa~DBl as 222''' <T< 226° and 
-3.4° < 6 < +0.7°. A host of features popul ate the 
CMa OBI association 



more than 30 nebulae (Herbst 
et al.|1978 ), the CMa Rl association containing a group 



ot stars ( mcluding HD 53623) embedded in a prominent 



reflection nebula (Van Den Bergh 



1966 



intercon nected H ii regions ot Sh 2-292, S 



I 



and the three 
2-296, and Sh2- 



297 see 



Grcgorio-Hcteml 2008 and references therein). 
Herbst S'Assousa (19771), based on the H i morphol- 



r (jweaver & Williams 1974), emission nebulosity, and 
1 dlst 



ogy 

the distribution & properties of high-mass OB sources in 
the region, postulated that a supernova remnant (SNR) 
could have induced star formation in CMa Rl. They es- 
timated the age of the supernova shell to be about half 
a million year s. The SNR hyp othesis was supported by 
the work of Herbst et al. ( 1978 1 - who estimated a similar 
age for most ot the stellar sources in CMa Rl us i ng op - 
tical and infrared studies; and that of Vrba et al. (1987), 
whose linear polarization survey of (jMa Rl is consis- 
tent with the model of supernova induced compression 
of an initially uniform magnetic field, though, it must 
be noted that stellar winds too can produce a shell like 
structu re. However, an alte r nate s tar formation hypoth- 
esis by [Reynolds fc Ogden ( 1978[ ) - based on the study 
of gas velocities and the tact that UV fluxes from two 
hot stars (HD 54662 and HD 53975) can account for the 
physical parameters of the shell - suggests that strong 
stellar winds or evolving H ii regions a re t he cause of 
star formation in t he region. |Blitz| ( |1980| ) and |Pyatuni1Ia| 
' &: Taraskin (|1986[) have s uggested the same 



Bonatto & Bica (20091 have investigated this region 
using I'wo Micron Ail Sky Survey (2MASS) data and 
concluded that it presents similar structura l properties 
as t ypical young and low mass open clusters. |Forbrich et| 



(|2004|) found a deeply embedded, very young source 
named as UYSOl (Unidentified 



intermediate mass. 
Young Stellar Object 1), that is driving a high-velocity 
bipolar CO outflow at the interface of the H ii region and 
the dark cloud. It was further resolved into two pr oto- 
stars, UYSOla and UYSOlb QForbrich et al]l2009l ). A 
recent Herschel PACS and SPIRf; mappmg ot the region 
around these sources has identified five cool and compact 
far infrared (FIR) sources, of masses of the order of a few 
Mf7), hidden in th e neighbouring dark cloud LDN1657A 
( |Linz et al.||2010| . 



Samal et al. 2007 2010 



egions dOjha et al. 


2004a|b|c 


Ojha et al. 20il|), 


m this pa- 



per, we explore the local star formation process around 
Sh2-297 region. The presence of multiple facets (like the 
ionizing star's radiation field, the nebulosity present in 




Fig. 1. — The optical color composite image of Sli2-297 region 
made using 6724 A [SII] image (red), 6563 A Ha image (green) 
and 5007 A [OIII] image (blue), observed using HCT (see Section 
|2A| . The black square marks the position of the ionizing star HD 

the region, and the dark cloud to the west) necessitates 
a multiwavelength study to decipher the star formation 
scenario. Hence, we carried out a detailed study of this 
region in optical, deep NIR, mid infrared (MIR) , and ra- 
dio wavelengths to learn about the central ionizing source 
and its effect on the region, the physical characteristics 
of the region, the spatial distribution of the YSOs and 
the nature of some selected individual YSOs. Based on 
the study of stellar population, their evolutionary sta- 
tus', and the morphology of the region, we have tried to 
deduce possible star formation activity. 

In Section[2j we describe our observations and the data 
reduction procedures. In Section [3j we have outlined the 
other data sets which were used in our study. We discuss 
the morphology of the region from radio, infrared (IR) 
and submillimetre data in Section [4j and the results per- 
taining to the stellar sources in Section [5] In Section [6) 
we discuss the star formation scenario in sh2-297 region. 
We present our main conclusions in Section [7] 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Optical •photometry 

Optical photometric observations were performed us- 
ing the 2 m Himalayan Chandra Telescope (HCT), 
Hanle (India), operated by the Indian Institute of Astro- 
physics, Bangalore. The telescope has Ritchey-Chretien 
optics with an f/9 Cassegrain focus. Both broad- 
band and narrowband observations were performed us- 
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ing Himalayan Faint Object Spectrograph and Camera 
(HFOSC) mounted on the telescope. HFOSC has a SITe 
(Scientific Imaging Technologies Inc.) 2k x 4k pixels 
CCD, with the central region of 2k x 2k being used for 
imaging. With a pixel scale of 0.3 , this translates to 
a field of view (FoV) of ~ 10 x 10 . Long and short 
exposure observations were taken in Bessell V (600 s 
& 20 s) and / (300 s & 10 s) filters (centered at q;2ooo 
« 07''05™13", (^2000 ~ -12°19'03") on 2009 October 21. 
[SII] (6724 A), Ha (6563 A), & [OIII] (5007 A) filter 
observations (centered at a2ooo ~ 07''05™17^, (52ooo ~ - 
12° 19 18 ), with an integration time of 300 s each, were 
done on 2010 January 18. Along with the object frames, 
several bias frames, twilight flat-fi eld frames, and the 
standard field SA98 ( Landolt] 1992 ) were also observed. 
The average seeing size ranged from ^2 to 3 full 
width half maxima (FWHM). The standard field was 
used to calculate the extinction coefficients so as to cal- 
ibrate our CCD systems. The astrometric calibration, 
with a position accuracy better than 0.3 , was done us- 
ing the 2MASS sources in the field. For all the filters, 
frames of equal exposure time were co-added to improve 
the signal-to-noise ratio. Photometry was carried out 
only for V and / filters. Image Reduction and Anal- 
ysis FacilitjF] (IRAF) data reduction package was used 
for the initial CCD image processing. Thereafter, us- 
ing DAOPHOT II package of MIDA^ the point spread 
function (PSF) photometry and the photometric mea- 
surements were performed. A variable PSF was made 
using several uncontaminated stars in the field. Finally, 
the magnit udes ( V and /) were calibrated in the manner 
outlined by Stetson ( 1987 ) , with a precision better than 
0.05 mag. In the present work, we have used only those 
sources which have magnitude errors < 0.15 mag (~ 3(t). 

2.2. Optical spectroscopy 

The optical spectroscopic observations were carried out 
using HCT and the 2 m lUCAA Girawali Observatory 
(ICO) telescope at Girawali near Pune (India), operated 
by Inter University Centre for Astronomy & Astrophysics 
(lUCAA), Pune. ' 

The HCT data was obtained on 2009 December 21 us- 
ing HFOSC, with the help of Grism 7 (3500 to 7000 A) 
- which has a resolving power of 1200 and a spectral 
dispersion of 1.45 A pixel^^. The spectrum of the star 

HD 53623 (02000 = 07''05™16.75^ (S2000 = -12°19'34.5"; 
V=7.99 mag) with an exposure time of 480 s, in addition 
to FeAr lamp arc spectrum (for wavelength calibration) 
and multiple bias frames, was obtain ed. The spectro pho- 
tometric standard star G91B2B (Massey et al.|1988|) was 
also observed with an exposure time of 1200 s. Spectral 
reduction was carried out using "APALL" task in IRAF 
reduction package. Final ly, th e flux-calibrated spectrum 



was obtained (see Figure [2(a) ) 

To cross-verify the weak diagnostic lines in the object 
spectrum from HCT, we also obtained a deeper spectrum 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy Inc. under cooperative agreement with 
National Science Foundation. 

^ MIDAS is developed and maintained by the European South- 
ern Observatory. 



(for 1200 s) for HD 53623 using ICO telescope, on 2009 
November 26. The telescope has a Cassegrain focus with 
a focal ratio of f/10. lUCAA Faint Object Spectrograph 
and Camera (IFOSC), with a 2k x 2k pixel CCD array, 
was used for the purpose. Grism 7, with a range of 3800- 
6840 A, resolving power of 1300, and a spectral dispersion 
of 1.49 A pixel" ^ was used. Multiple bias frames, halo- 
gen flat-field frames, and HeNe lamp arc spectrum (for 
wavelength calibration) were obtained along with the ob- 
ject spectrum. The reduction was done using "DOSLIT" 
task in IRAF. No flux calibration was done for this case. 
Fin ally, t he normalized spectrum was obtained (see Fig- 
ur ej2(b)| ). 

lb identify the Ha emission line stars, two grism slit- 
less spectra (420 s each) - using Grism 5 (5200-10300 
A) and Ha broadband flher (A = 6563 A, AA = 500 
A) - were also obtained with HCT on 2007 November 16 

(centered at a2ooo ~ 07''05™23^ ^2000 ~ -12°19'00"). As 
a reference, for the identification of stars whose slitless 
spectra were obtained, two Ha broadband filter images 
of the field were taken as well. Here, the first step was to 
average the two photometric Ha filter images and the two 
slitless spectra of the region to improve their signal-to- 
noise ratio. The world coordinate system (wcs) coordi- 
nates were then implemented on the averaged photomet- 
ric image by identifying a bright set of stars in the FoV, 
from the USNO catalogue, and finally synthesizing their 
J2000 coordinates and pixel coordinates from the im- 
age, using the generic IRAF commands "CCMAP" and 
"CCSETWCS". The averaged spectrum was then ana- 
lyzed in ds9 image display device, and the stars with an 
Ha emission line (detectable as an enhancement over the 
continuum) were identified manually by careful visual in- 
spection, with the help of photometric image of the field. 
19 Ha emission line stars were identified. The detection 
limit of this identification is ^ 3.3 A in terms of equiva- 
lent width, with the faintest star having a V magnitude 
of 20.07. Fifteen were located in our NIR FoV, and hence 
the NIR magnitudes for the remaining 4 were taken from 
2MASS catalogue for further analysis. 

2.3. Near infrared observations 

The region of interest, being a star forming region, has 
considerable nebulosity and dark cloud to the west, and 
hence is obscured at optical wavelengths. Moreover, we 
need to identify and classify the embedded YSOs which 
emit at longer wavelengths, hence we use NIR observa- 
tions. NIR observations, centered on a2ooo ~ 07''05™11'', 
(52000 « -12°18'5l" for the object, in J (A = 1.25 ^im), H 
(A = 1.63 /im) and (A = 2.14 ^m) bands, were carried 
out on 2007 February 27 using the f/10 Cassegrain 1.4 m 
Infrared Survey Facility (IRSF) telescope. South Africa 
with SIRIUS (Simultaneous InfraRed Imager for Unbi- 
ased Survey) - a three color simultaneous camera. The 
camera is equipped with three Ik x Ik HgCdTe arrays, 
with a FoV of ^ 7.8 x 7.8 , and a pixel scale of 0.45 . 
Further information about the instrument can be found 

For each 



and Nagayama (2003) 



in jNagashima ( ,1999) _ ^ 

band, 4U tranies eacn, centered at 9 dithered positions, 
were taken. The integration time for each frame was 10 
s, giving a total integration time of 3600 s (9x40x10) in 
each band for the Sh2-297 region. The sky condition was 
photometric, with the average seeing size ranging from 
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Fig. 2. — (a) Flux-calibrated spectrum of the central ionizing 
star, HD 53623, obtained with HCT (b) Normalized spectrum of 
HD 53623 obtained with IGO. 

1.35 to 1.60 FWHM during the observations. 

The standard data reduction procedure - involving bad 
pixel masking, dark subtraction, flat field correction, sky 
subtraction, combination of dithered frames, adding wcs 
coordinates to the image and PSF photometry - was car- 
ried out using IRAF. Astrometric calibration was imple- 
mented using the 2MASS sources in the field. About 
30-40 stars, with good profiles, were selected for calibra- 
tion for each band, and an accuracy of better than ± 
0.04 was achieved for our absolute position calibration. 
After the initial processing, photometry was carried out 
on a final FoV of - 7.5' x 7.5'. The PSF photometry 
was done using "ALLSTAR" algorithm of "DAOPHOT" 
package in IRAF. Around 11 to 13 isolated bright stars 
were taken to make a PSF for each band. Subsequently, 
the final magnitudes were calibrated using the 2MASS 
stars in the field, with an rms of ~ 0.07 mag for all the fil- 
ters. To ensure good photometric accuracy, we restricted 



our catalogue to sources having uncertainty < 0.15 mag 
(~ 2(t), independent of brightness, in all the three bands. 
Also, it was found that the sources with Kg magnitude 
< 8 were saturated in our catalogue; hence their mag- 
nitudes, in all three filters, were replaced by the corre- 
sponding 2MASS values. We evaluated the completeness 
limits for our images by doing artificial star experiments. 
Stars of varying magnitudes were added to the images, 
and then it was determined what fraction of stars were 
recovered in each magnitude bin. The recovery rate was 
greater than 90% for magnitudes brighter than 17, 17, 
and 16 in J, H, and Kg, bands, respectively. The re- 
covery percentage fell to 85% for Kg magnitudes in 16 
to 16.5 mag bin. The change in completeness limit is 
negligible on using the photometric error cutoff at 0.15 
mag. 

The final calibrated catalogue was analyzed using the 
CC and CM diagrams (see Section [5]). For this pur- 
pose, the required magnitudes and colors were trans- 
formed into CIT (California Institute of Technology) sys- 
tem using the color transformation equations given on 
the CALTECH websitc[5 In addition to the target re- 
gion, a control field centered at a2ooo ~ 07''06™20^, (520oo 
w -12°01 40 (about 24 to the north-east of our object 
field) was also observed. It was reduced in a similar man- 
ner as the object frames. The control field was used to 
assess the level of contamination from foreground sources 
in the target field. It was used to delineate the maximum 
H — K color for the foreground sources, so as to iden- 
tify the young stellar sources in the region - which will 
have a larger color by virtue of extinction and the pres- 
ence of circumstellar disk which results in excess K-hanA 
emission. Figure [3] shows a false NIR color composite (J : 
blue, H : green, iTs : red) image of Sh2-297 region, which 
reveals nebulosity around the central and south-eastern 
region, and the dark cloud towards the west. 

2.4. Radio continuum observations 

1280 MHz observation was carried out on 2007 July 
12, while that for 610 MHz on 2007 July 19, using Giant 
Metrewave Radio Telescope (GMRT) array. The GMRT 
array consists of 30 antennae, arranged in an approxi- 
mately Y-shaped configuration, with each dish being 45 
m in diameter. 12 antennae are located randomly in a 
compact 1 kmxl km central square area, while the re- 
maining 18 are located along the three radial arms (6 
along each arm), each arm length being ~ 14 km. The 
maximum baseline is ~ 25 km. Further details can be 



found in Swarup et al. ( 1991 ). For our observations, VLA 
phase cafibrator "07^^-175" and flux calibrator "3C147" 
were used, at both the frequencies. The details of the 
observations are given in Table [l] 

The radio data from GMRT was analyzed with the 
help of Astronomical Image Processing Software (AIPS) 
distributed by NRAO. The data was edited to flag the 
bad antenna pairs, baselines, and time ranges using the 
"VPLOT" , "UVFLG" , and "TVFLG" tasks. The final 
image was produced using a set of tasks which involved 
fourier inversion, cleaning (using the task "IMAGR" in 
AIPS), followed by a few iterations of phase-only self- 
calibration (till the rms convergence) to remove the ef- 

http:/ /www. astro. caltech.edu/~jmc/2mass/v3/transformations/ 
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TABLE 1 
Details of GMRT Observations 



1280 MHz 



610 MHz 



Date of Obs. 
Phase Center 

Flux Calibrator 
Phase Calibrator 
Cont. Bandwidth 

Primary Beam 
Resolution of maps 

used for fitting 
Peak Flux Density 
rms noise 
Integrated Flux Density 



2007 July 12 

aaooo = 07'^05™ 12.30" 
<S2000 = -12°19' 01.14" 
3C147 
0735-175 
16 MHz 
24' 

23.5" X 13.9" 
19.42 mjy/beam 
2.32 mJy/beam 
0.67 Jy 



2007 July 19 
02000 = 07''05™12.30= 
52000 = -12°19'01.14" 
3C147 
0735-175 
16 MHz 

43' 

23.8" X 13.6" 
41.33 mJy/beam 
3.24 mJy/beam 
0.64 Jy 



.25-0 20.0 . 15.0 10.0 05.0 7:-05:()0.0 



3C147). The temperature for the source will be T£ 



RA (J2000) 



Fig. 3. — Color composite made using J (blue), H (green) and 
Kb (red) images. White square marks the position of ionizing star 
HD 53623, white cross shows the position of IRAS 0702 9-1215. and 
the mage nta plus marks the location of UYSOl from [Fbrbrich et| 

p^lpoolt . 

fects of the atmospheric and ionospheric phase distor- 
tions. 

Also, since our source is located near the Galactic 
plane, we had to take into account the system tempera- 
ture corrections for GMRT. Since, at metre wavelengths, 
there is a large amount of radiation from the Galac- 
tic plane, the effective temperature of the antennae in- 
creases, which has to be corrected for. Thi s was done 
by obt aining the sky temperature map of |Haslam et al.] 
( 1982 1 at 408 MHz, and rescaling the image by a scal- 
ing factor which is just the ratio of system temperature 
towards the target source and the flux calibrator (here 



T; 



sys 



where T 



T408 x(f/408)~ ^-^ (using the spectral 
119821), with T408 being the 



Haslam et al. 



index given in 

temperature value from Haslam's map at 408 MHz and 
"f" being the frequency (in MHz) for which the correc- 
tion is to be calculated. Tgys, the system temperature, 
was obtained from the National Centre for Radio Astro- 
physics (NCRA) website]^ for both the frequencies. The 
temperature for the flux calibrator will merely be Tgys, 
hence the scahng factor is (T^ -I- Tgys ) /Tgys- 

In addition, for 1420 MHz, the archival image from 
Very Large Array (VLA), for the observation date 1998 
December 14, was also used (Project ID AF346). 

Radio continuum interferometric data, at 610 MHz, 
1280 MHz and 1420 MHz was used to probe the ionized 
gas component. The GMRT maps were convolved to 
a resolution of about 24 xl4 , approximately same as 
the VLA map. After obtaining the contour maps for all 
the three frequencies, the integrated flux was found by 
fitting a 3cr contour level box and using the AIPS task 
"IMSTAT" (see Table [l]). Using these radio continuum 
data points, the various phys ical parameters of the region 
were calculated (see Section 4.2). 



3. OTHER AVAILABLE DATA SETS 

3.1. Near infrared data from 2MASS 

We used the 2MAS^ data in the field for our photo- 
metric and astrometric calibrations. The stars were ob- 
tained using the task "IMTMC" in "WCSTOOLS" pack- 
age in IRAF, for the J (A = 1.25 ^im), H {X = 1.63 ^im) 
and Ks (A = 2.14 /xm) bands. All the sources which were 
obtained for calibration had a good quality flag value 
of "AAA" for the magnitudes. We also substituted the 
2MASS magnitudes for those sources which were found 
to be saturated in our IRSF catalogue, i.e. those for 
which we had IRSF Ks magnitude < 8. 

3.2. Near and mid infrared data from WISE 

Wide-field Infrared Survey Explorer (WISE) surveyed 
the entire sky in the NIR bands at 3.4 /im (wl) and 4.6 
fim (w2), and in the MIR bands at 12 //m (w3) and 22 

* http:/ /www.ncra.tifr.res.in/~gtac/GMRT-specs.pdf 
^ This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Center, 
California Institute of Technology, funded by NASA and the NSF. 
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fim ( w4), with an astrom etric precision better than ~ 
0.15 (Wright et al. 12010 1. The spatial resolution was ^ 

II II 

6 for the first three bands and ^ 12 for 22 /im. 

Since the WISE (Preliminary Release April 2011) cat- 
alogue had many sources missing (e.g. for wl band, only 
~ 35% of the sources detected in our photometry were 
present in the release catalogue), as could also be visually 
discerned, we carried out aperture photometry on the 
archival WISE images - for the same FoV as o ur IRSF 
NIR field - following the method described by |Koenig| 



et al.| ( |2012[ ). This significantly improved the number 
ot sources detected. Photometry was done using the 
DAOPHOT package in IRAF; with the aperture radius, 
inner radius of the sky annulus, & outer radius of the 
sky annulus being 5 ,7.5 , & 25 for the 3.4 & 4.6 ^m 
bands, and 7.5 ,10 , & 25 for the 12 /im band, respec- 
tively. The 22 /im band was not used mainly because of 
its low resolution, as well as the dominant nebulosity and 
the presence of not many point sources (12 sources were 
detected out of which 9 were in highly nebulous central 
region). The photometry for each band was calibrated 
with the WISE preliminary release catalogue. Care was 
taken not to use the objects flagged as "D" (diffraction 
spikes) , "H" (halos from bright sources) , "P" (persistent 
artifacts), or "O" (optical ghosts of bright sources) in 
the release catalogue for calibration. For the final WISE 
catalogue, we took only those sources which had all mag- 
nitude errors < 0.3, to make sure that they had reliable 
photometry. 

Our IRSF NIR catalogue was subsequently matched 
with the WISE catalogues within a 3 matching radius. 
187, 124, and 9 sources detected in our NIR FoV had 
3.4 /im, 4.6 /im, and 12 /im counterparts, respectively. 
The resultant final catalogue was used to get the SEDs 
of various sources, and hence in constraining the mani- 
fold physical parameters l ike a ge, mass, accretion rates 
etc. for them (see Section 5.6 1. The 3.4 /tm and 12 /tm 
bands contain prominent poly cyclic aromatic hydrocar- 
bo n (PAH) features at 3.3, 11.3 & 12.7 /im ( [Wright et 
al.||2010| iSamal et al. 112007') in addition to t he contin- 



uum emission, and hence can be used to get an idea of 
the photodissociation region. The 22 /tm band can be 
used to examine the warm dust emission - the stochastic 
emission from small gr ains as well as the thermal emis- 
sion from large grains ( [Wright et aLl[2010l ). 



4. MORPHOLOGY OF THE REGION 

4.1. The ionizing source 

The ionizing source of Sh2-297 region is the star HD 
53623, associated with the reflection nebula CMa Rl 



^an Den Bcrgh'1966|). |Claria| ( |1974[ ) and |Houk fc Smith- 
[Moorc (1988) estimated the spectral type ot this source 
as~BTVnand BlII/III, respectively, using objective- 
prism plates from Michiga n Curtis-Schmidt telescope, 
CTIO. Herbst et al.| (1978) have proposed the spectral 
type as BO. 5 IV-V using spectrogram analysis. Here we 
use the optical spectroscopic data to find out the spec- 
tral type of this source. The reduced spectra from HCT 
and IGO, for the wavelength range 3950 to 4750 A, have 
been shown in Figure [2] The spectra w ere analyzed with 



A, and low strength at 4686 A suggest that the star does 
not belong to 0-type spectral class. The most prominent 
lines in our spectra are the He i lines at 4009, 4026, 4144, 
4387, 4471, and 4713 A. The Balmer lines at 3970 A 
(He), 4102 A (H^), and 4341 A (H7) are clearly seen too. 
We also see C iii -I- O 11 blend at 4070 and 4650 A. Weak 
Si IV line at 4116 A, weak O 11 lines at 4415-4417 A, and 
absence of Si ill lines indicate that the star most probably 
does not belong to giant luminosity classes. The prepon- 
derance of neutral helium lines and the comparison with 
the standard Walborn catalogue suggest a most likely 
spectral t ype of B0V/B0.5V for the star. |Tjin A Djie et 
al. (2001) have suggested, using Ha emission study, that 
tms source does not contain any circumstellar disk, most 
probably due to photoevaporation by UV photons or due 
to influence of some supernova in the nearby region. 

4.2. Physical properties of the region 

High resolution studies are important as they help us 
in understanding the structure of a region. Sh2-297 has 
been observed with very coarse resolutions in previous 
radio surveys, which makes it harder to recognize the 
subtler feat ures in the radio emissio n. For example, the 

IT 



survey by Felli & Churchwell 
cally identified H 11 regions at 



(h972) observed 168 opti- 
raODMHz, including Sh2- 
297, with a spatial re solution of 1 . Similar observations 



were carried out by Pyatunina (1980) for wavelengths 
ranging from 2 to 11 cm, with resolutions of ~ 1 —4.3 . 
At larger fre quencies (10 GHz), th e CMa Rl region was 
obse rved b y Nak ano et al. (1984) with a resolution of 

2.7' 



Fich 



(11993), in a VLA survey, observed the region 
at 1.4 GHz with 43" resolution. Figure [i] shows GMRT 
high resolution greyscale contour maps, for 610 MHz (^ 
10" X 5") and 1280 MHz (- 7"x 6"), of Sh2-297 H 11 re- 
gion. The intricate morphology of Sh2-297 - showing the 
complex structures and the clumpy nature - can very well 
be discerned on examination. The radio emission mor- 
phology from our observation as well as the location of 
H II region at the tip of the molecular cloud (see Section 
4.3 1 suggest that the ionized gas is possibly undergoing 
a champagne flow (also known as blister region) ( |Teno-| 
rio Tagle et al.||1979l lAUoin & Tenorio Taglel 19791), with 



the extended emission towards the east and the head of 
the ionized champagne flow located to the west towards 
LDN1657A (see Section |6|). 

The radio continuum observations of the region using 
GMRT, and archival data from VLA, were used to cal- 
culate the various physical parameters which character- 
ize an H II region (emission measure, electron density, 
stromgren radius, and dynamical timescale). Assuming a 
spherically sy mmetric and homogen e ous i onized region, 
the model of Mezger fc Henderson| ( |1967) for free-free 
emission was used. According to this, the flux den- 
sity arising du e to the region is g iven by the following 
(adapted from Mezger et al.||1967 Equations 1 and 3) : 



5"^ = 3.07 X lQ-'^T^v'^VL{l - e^^^")) 



(1) 



the help of W alborn spectra catalogue (Walborn & Fitz 



patrick|1990). The absence of He 11 lines at 4200 & 4541 



t{v) = 1.643 X Iff'aT-^-^'^v-^-^nll (2) 

where. Si, is the integrated flux density in Jansky (Jy), 
Te is the electron temperature of the ionized core in K, 
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RIGHT ASCENSION (J2000) 



100 200 500 1000 2000 5000 

Frequency ( MHz) 

Fig. 5. — Observed and fitted flux density (solid line) for the 
H 11 region. Filled circles and triangle represent GMRT & VLA 
data, respectively. The flux densities were calculated for an ap- 
proximately same resolution of ~ 24 X 14 . 




RIGHT ASCENSION (J2000) 



Fig. 4. — (upper panel) High resolution GMRT radio map for 610 

MHz {~ 10"x5"). Contours have been drawn at 12, 16, 20, 24, 
& 30 % of the peak flux (28.30 mjy/beam). (lower panel) High 

resolution GMRT map for 1280 MHz (~ l" X g"). Contours have 
been drawn at 10, 15, 20, 25, 30, 40, & 50 % of the peak flux (10.47 
mJy/beam). 

V is the frequency in MHz, Ue is the electron density in 
cm"^, / is the extent of the ionized region in pc, r is 
the optical depth, a is the correction factor, and Q is 
the solid angle subtended by the beam in steradian (ap- 
proximating the beam with a Gaussian, the beam solid 
angle is given by = 1.133 x 9rnaj x 6',„„i, where 9maj 
and 9min are major and minor half power beam widths). 
The factor ri^l denotes the emission measure (in cm^^ 
pc), a measure of optical depth of the medium. For our 
calculations, we set the valu e of a as 0.99 (obtained from 
Mezger fc Henderson|ri967[ Table 6), and th at of T^, to 
lUUUU K - apphcable to classical H ii regions (Stabler & 



Palla 2004). Thereafter, treating the emission measure 
as a tree parameter and using the three data points (flux 
density vs. frequency) from the three maps, a non-linear 



regres sion was implem ented, following a procedure simi- 
lar to Vig et al. ( 2006 ) (see Figure [S]) , which yielded the 
value ot the tree parameter (emission measure), for the 
best fit, as 9.15 ± 0.56 x 10^ cm^'' pc. The extent of 
this extended H ii region is ^ 5 , as obtained from the 
VLA image. At a distance of 1.1 kpc, this translates to 
an extent of 1.6 pc. From the emission measure and the 
extent, the electron density is obtained to be 756 ± 46 
cm~^. The low values of the emission measure and the 
electron density reinforce the fact that this is a classical , 
and thus, evolved H ii region (see Stahler fc Palla|2004 ). 

Further, assuming a homogeneous and spherically sym- 
metric nature for the H ii region, the luminosity of the 
Lyman continuum photons (in ph otons s~^) fro m the 
central ionizing source is given by ( Moran 1983 , Equa- 
tion 5) : 



5* = 8x10 



43 



Si/ 
mjy 



V 104/^ / 



-0.45 



(- 

\ kpc 



GHi 



0.1 



(3) 

where Si, is the integrated flux density in mJy from 
the contour map, D is the distance in kiloparsec, Tg is 
the electron temperature, and v is the frequency in GHz 
for which the luminosity is to be calculated. For our cal- 
culations, we used the 1420 MHz flux value (0.76 Jy). D 
was taken to be 1.1 kpc and Tg as 10'' K. The value for 
S** was determin ed to be 7.62x10 ' ^^ pho tons s~^ (hence 



log 5* is 46.88). Bonatto & Bica (2009) have estimated 



the main sequence (Mb) age of HU 53623 as 5 ± 3 Myr 
therefore assuming a l uminosity class V for it, a com- 
parison of log 5* with Panagia' ("1973, Table H) yields 
the spectral type as BO, in agreement with our optical 
spectroscopic analysis. The lower value of our calculated 
luminosity ( log Si, ) as comp ared to the theoretical value 



(47.63) from Panagia ( 1973 ) could be because of absorp- 
tion of photons by dust m the region. 

As the central massive star ionizes the ambient medium 
around it, the ionization front moves outwards, and the 
whole H II region expands till an equilibrium is reached 
between the number of ionizations and recombinations. 
The underlying assumption here is that the ambient 
medium has homogeneous density and temperature. Us- 
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ing the formulation given by Stromgren ( |1939 ), the ra- 
dius of the H II region at this point - called Stromgren 
radius - is given by : 



Rs = 



yman 



V 47^2/32 



1/3 



(4) 



where Rg is the Stromgren radius (in cm), Uo is the 
initial ambient density (in cm~'^), and (32 is the total 
recombination coefficient to the first excited state of hy- 
drogen. The value of /32 , for a of 10"^ K. is taken 
to be 2.6x10-13 cm^ s'^ ( [Stahler fc Palla||2004| ). We 
est imated the initia l ambient density as tollows (similar 



to [Stromgren 1939 ) . From the luminosity calculations 
for Lyman continuum photons, we have obtained a cal- 
culated value of 7.62x10*^ photons s^^ (5'*), which is 
about 18 % of the theor etical output of 4.27x10'*'' pho- 
tons s"* ( Panagia||l973 ). This implies that rig — O.lSrio, 
hence yielding Uo as 4202 ± 253 cm-"^, using the above 
calculated value of the electron density. Using these val- 
ues for 5* and Uq, the stromgren radius turns out to be 
0.051 ± 0.002 pc. 

As the gas expands, it sends a shock front into the 
neutral medium. This shock front precedes the ionization 
fro nt. At this sta ge, the radius of the H ii region is given 
by ( |Spitzer|[l978l ): 



(5) 



where R{t) is the extent of the H ii region at time t, 
and cjj is the speed of sound, taken to be 11x10^ cm 
(Stahler & Palla 2004 1. Using the current extent of 1.6 
pc as the value ot R(t), and the above calculated value 
of Rs in this equation, we obtain the value of t, referred 
to as the dynamical timescale, as 1.07 ± 0.31 Myr. 

It should be noted that, to calculate Rg and t, we have 
estimated Ho (4202 ± 253 cm'^) assuming a luminosity 
class V for the massive star, and that only part of the 
radio flux is used for ionization. Also, we have used the 
current value of electron density. However, at the begin- 
ning of massive star formati on. Un could be of the or der of 



lO^'-lO^ cm-3 ( (Ward-Thompson fc Whitworth|201l[) . Be 



sides, the expansion velocity decreases with time. Mence, 
the electron density and the dynamical timescale calcu- 
lated here should be treated as conservative lower limits. 

4.3. Infrared and submillimetre structures 

Figure § shows the i2cO(J=3-2) image of the region 
(public processed data for observation date 2011 Febru- 
ary 20, Project ID: MllAH48Af] from James Clerk 
Maxwell Telescope (JCMT). TheHlSX 8.28 /xm con- 
tours (in yellow), SC UBA 850 /xm contours (in red) (Di 
Francesco et al.|2008 ), and the 1280 MHz radio contours 
from GMRT (in white) are overlaid. The MSJ<QA band 

® The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre on behalf of the Science and Technol- 
ogy Facilities Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research 
Council of Canada. The archival data was downloaded from 
http: / /www. cade. hia.nrc.gc.ca/jcmt / search/product / 

' I'his research made use of data products from Midcourse Space 
Experiment. Processing of the data was funded by Ballistic Mis- 




RA (J2000) 



Fig. 6. — The ^^CO(J=3-2) image (log greyscale with intensity 
range -1 to 30 K km s-l) of the region from JCMT, overlaid with 
MSX 8.28 ^im contours (yellow) showing the PAH distribution, 
SCUBA 850 /im contours (red) which trace the cold dust, and 
1280 MHz radio continuum contours (white) indicating the ionized 
gas. The CO image reveals a maximum column density of 87.28 
K km s~^. MSX and SCUBA contours have been plotted at 10, 
15, 20, 25, 30, 40 & 50 % of their peak values, which are 49.666 
mW/m'^-Sr and 1.619 Jy/beam, respectively. Radio contours (~ 

25 resolution) have been plotted at 3, 6, 9, 15, and 21cr (cr denotes 
the rms noise, which is 1.03 mjy/beam) levels. The blue square 
marks the ionizing star HD 53623, blue plus symbol marks the 
position of UYSOl from Forbrich ct al. (2004), blue circle marks 
the position of C2H peak from Beuthcr et al., ^2008,). blue cross 



denotes IRAS 07029-1215, and the sources from 
have been labelled and marked in white circles. 



rinz et al.| | |2010[ l 



(A = 8.28 /im; containing PAH features at 7.7 & 8.7 /.tm 
besides continuum emission) contours give an idea of the 
P AH dis tribution in the region (see e.g. Deharveng et 



alj2005|. i2co(J=3-2) traces the molecular hydrogen 
while SCUBA 850 ^m traces the cold dust. As can be 
seen from the figure, there is a clear correlation between 
the ionized gas and 850 /im contours in Sh2-297 central 
region. The ^^CO image shows a D-shaped ring to the 
south of HD 53623, which is also traced by the MSX A 
band contour s, and is faintly visibl e in narrowband H2 
S(l) image of Forbrich et al. (2009 see their Figure 1). 
This is probably due to the sweeping up and piling of 
material from the expanding H 11 region, and as such it 
is a potential site for future triggered star formation by 
coUect-and-coUapse process. Also noticeable is the fact 
that, towards the west, PAH extends roughly upto the 
interface of the H 11 region and LDN1657A. 

sile Defense Organization with additional support from the NASA 
Office of Space Science. 
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The WISE w4 band which traces the warm dust and 
thermal emission from protostars (not shown here), and 
SCUBA data show that the warm and cold dust extend 
well beyond the H ii region, and also towards the cold, 
dark cloud LDN1657A. The gas and dust correlation is 
evident in the FIR IR AS 100 /im wavelength as well 
( Gregorio-Hetem|2008 ) . The nebula around the star HD 
53(j 23 sho ws a high level of FIR emission. |Beuther et 
ah] ( p08| found a C2H peak near IRAS 07029-1215, co- 
mciding svith UYSOl (see Figure |6|, t hereby iiidicating 
that this is a young star forming region ( Vasyunina et al 
2011 1. |Kim et al.| ( |2004| ) have probed the dense molec- 
ular gas using the optically thin J = 1 — transition of 
One of the peaks is at ^ = 225.47°, b = -2.80" 
(outside the range of Figure |6]), showing the presence of 
dense gas further towards the south-west of Sh2-297 as 
well. 



5. STELLAR POPULATION IN THE REGION 

5.1. NIR color-color diagram 

The NIR CC diagram for the IRSF sources, sh owing 
the J — H vs. H — K colors, is given in Figure 7(a) 



The thick red line denotes the locus of MS stars, the ma- 
genta line shows the locus of giant stars while the dot- 
ted black straight line is the locus of Classical T Tauri 
Stars (CTTS). The loc i for the MS and giants were taken 
from [Bessell fc BrettI f l988) while those for the CTTS 
from Meyer et al. 1TI397). The three parallel dashed and 
slanted lines are the reddening vectors [Aj/Ay = 0.265, 
Ah/Av = 0.155 and Ak/Av = 0.090) for the CIT sys- 
tem from Cohen et al. (1981). The crosses mark every 
5 mag increase in Ay. We converted the colors from 
2MASS to CIT system, as was discussed in Section [273} 
We used the CC diagram to identify and study the na- 
ture of YSOs in the region. Since young sources have 
a significant infrared excess, which in turn can be used 
to identify them, we divided our CC diagram into thre e 
regions - F, T and P (similar to Ojha et al. 2004b|c 



The region "F" mostly contains the held stars and Class 
III sources (weak-li ne T Tauri stars). Th e region "T" 
contains the CTTS ( |Lada fc Adams||1992| ), Class II ob- 
jects with large NIR excesses, and reddened early type 
ZAMS (Zero Age Main Sequence) stars with ii'-band ex- 
cess emission. The "P" zone sources are most likely Class 
I objects (protostar-like objects) with circumstellar en- 
velopes. The NIR CC diagram is likely to be contami- 
nated by the foreground and background field stars. To 
assess this, we als o plot ted the CC diagram for the con- 
trol field (Figure 7(b)). On comparison, we find that 
almost all our control field sources are confined to the 
left of the middle reddening vector, and a few which are 
to its right are below the CTT locus. Thus, the "T" 
and "P" regions of the object field are mostly uncontam- 
inated by the field stars, and hence the sources in these 
regions can be assumed to be YSOs. 

As there is considerable nebulosity in the region, it 
must be kept in mind that there could be dust knots 
emitting at near and mid infrared wavelengths. These 
could be mistakenly identified as Class I/II sources in the 
CC diagram. Hence, this CC diagram has its limitations 
when it comes to the identification of the YSOs, and 
the information from it is only complementary to other 
techniques like spectroscopy, to understand the nature of 
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Fig. 7. — (a) NIR CC diagram for the object field. The red curve 
shows the MS locus and the magenta curve shows the locus of giant 
stars. The locus of CTTS has been shown by a black dotted line. 
The three parallel and slanted lines are the reddening vectors, with 
crosses on them marking an interval of Av= 5 mag. Blue circles 
denote the Class Il-type sources, orange crosses the Ha emission 
line stars, and the green square the central ionizing star, (b) The 
CC diagram for the control field. 
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objects. The spatial distribution of the YSOs identified 
was used to pin down the st ar fo rmation process going 
on in this region (see Section 5.7). 



5.2. NIR color-magnitude diagram 

The NIR CM diagrams, K/{H - K), of the region 
and the control field are shown in Figure Is] The nearly 
vertical solid lines are the loci of ZAMSreddened by 
Ay = 0, 15, 30 and 45 mag. The slanting lines are the 
reddening vectors for the corresponding spectral types. 
A comparison of the object field and the control field 
CM diagrams shows an apparent separation between the 
stellar sources at around H — K ^ 0.6. Our control field 
shows almost all sources having H — K < 0.6, thereby 
indicating that sources with H ~- K < 0.6 are most likely 
the foreground stars. The sources lying to the right of 
the H — K ^ 0.6 c ut-off line (shown with a dashed ver- 
tical line in Figure 8(a) ) have an excess emission in the 
i^-band and can thus reasonably be assumed to b e young 
protostars. A CC diagram {Ks — 3.4/3.4 — 4.6, 'Koenig 



et al. 2012) for the red sources which had WISE coun 
terparts (not shown here) showed their location in the 
diagram to be consistent with the location of Class I 
sources. The central ionizing star has been marked with 
a green square on this diagram. The NIR CM diagram 
indicates a Bl spectral type for the ionizing source, which 
is within one sub-class of that estimated from the optical 
spectroscopy and radio continuum observations. How- 
ever, a spectral type derived from just two data points 
{H and K) is not very reliable because of uncertainty in 
distance determination. It just serves to give an idea of 
the spectral type, the spectrum analysis being canonical 
in this respect. 

We estimated the extragalactic contamination (ex- 
tincted or faint background galaxies) in our sample of 
NIR excess sources. We assume that the extragalactic 
sources can be seen through maximum extinction, which 
for our case is Ay ~ 25 mag (A^ ^ 2.25 mag) (see 
Section 5.3). Accounting for the interstellar extinction 
and converting the i^g magnitude to AB system (since 
K p, counts are tabu lated in AB photometric system in 



Keenan et al. (2010)), the total number of galaxy counts 
in our NIR region, i.e. ~ 56.25 arcmin^, turns out to be 
6 ± 2 down to our faintest detected magnitude of ^ 17 
in the Ka, band. 

A list of all the YSOs, with their respective magnitudes 
in optical and IR, has been given in Table [2j 97 YSOs 
were analyzed in total, with 93 of them in our NIR FoV 
of ~ 7.5 X 7.5 , and 4 Ha emission line stars outside this 
FoV (but considered for analysis nevertheless). Out of 
the 93 YSOs, 15 are Ha emission line stars detected using 
grism slitless spectroscopy with a V limiting magnitude 
of 20.07, 76 are NIR YSOs (26 Class Il-type sources + 50 
red sources with H~K > 0.6) with a limiting magni- 
tude of - 17, and 2 are LDN16 57A-2 and LDN1657A-3 

2010[ ) (only LDN1657A-3 was 
-like and H - 



sources from Linz et al. 
detected in NIR). Class 



K > Q.Q sources 
which were also detected as Ha emission line sources were 
classified as Ha emission line sources. The H — K > 0.6 
sources which overlapped with Class Il-like sources were 
classified as Class Il-like sources. 

5.3. Extinction in the region 
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Fig. 8. — (a) NIR CM diagram for our object field. The nearly 
vertical solid lines are the loci of ZAMS stars reddened by Av=0, 
15, 30, & 45 mag. The parallel, slanting lines denote the reddening 
vectors. Blue circles denote the Class Il-like objects, orange crosses 
the Ha emission line stars, and the green square the central ionizing 
star. The H — K cut-off at 0.6 has been marked with a dashed 
vertical line, (b) CM diagram for the control field. 
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TABLE 2 

Optical & infrared magnitudes for the 
YSOs 

i'able 2 has been given at the end ot the paper 



The average extinction in the region was calculated 
using the CC diagram as well as the extinction map of 
the region. 

To calculate the average reddening using the CC dia- 
gram, we used the sources in the "F" and "T" zones of 
our diagram. The sources in the "T" region were dered- 
dened to the CTT locus, while the sources in the "F" 
region were dereddened to the MS locus which was ap- 
proximated by a straight line asymptote. This gives the 
color-e xcess for each sourc e. Subsequently, the reddening 
laws of Cohen et al. ( 1981 ) were followed to calculate the 
Ay for each source. Finally, we fit the set of Ay values 
(ranging from ^ 0.02 to 16 mag) obtained for each of the 
sources, with a normal distribution. The mean of this 
distribution was found to be 2.70 mag. 

The average visual extinction for the region was also 
found by making an extinction map ("NI CE" method, 
Lada et al.||1994[ [Kainulainen et al.||2007| . To do this, 
we hrst obtained the intrinsic H — K color from the con- 
trol field. Thereafter, the color-excess, E(i? — K), and 
then Ay for each of the sources in our catalog was cal- 
culated. The color-excess is related to visual extinction 
by : E(i? — K) = 0.065x Ay, which can be calculated 
from the reddening laws. Ay values ranged from 0—25 
mag, as we took all sources with at least H and K mag- 
nitudes, and thereby even redder ones (undetected in J), 
here. Subsequently, the extinction map was made, where 
the extinction at each pixel was calculated by a weighted 
Gaussian mean around it. For our case, we took a Gaus- 
sian with an FWHM of 1 . Finally, the mean pixel value 
for our map was calculated using the task "IMSTAT" in 
IRAF, which gave us an average visual extinction of 2.69 
mag. 

Hence, we find similar Ay values from both the meth- 
ods. We used 2.70 mag as our value for Ay in the optical 
CM diagram to e stimate the mass and age of the YSOs 
(see Section 5.4). It gives the mean extinction of the 
field. However, owing to the nebulosity in the central 
region and the presence of dark cloud to the west (see 
Section 5.7), individual extinction for the sources will 



vary depending on their location. 



5.4. 



Optical color-magnitude diagram 

Figure |9] shows the V/V-I CM diagram for the YSOs 
in the region. Here, Ha emission line sources have been 
shown in orange crosses. Class Il-like sources in blue cir- 
cles, and sources with H—K > 0.6 have been represented 
by red circles. Most of the sources detected at optical 
wavelengths are Ha emission line sources, with only three 
Class Il-like sources and one source with H — K > 0.6. 
This is expected as Ha emission line sources are the most 
evolved YSOs in the field. I n Figure[9} the isoch rone for 2 
Myr for solar metallicity by Girardi etral.| (2002), and pre- 
main sequence (P MS) isocnrones tor 0.1 to 5 Myr from 
Siess et al. (20001 have been overlaid. The isochrones 
were corrected tor dista nce ( 1.1 kpc) and reddening (Ay 




Fig. 9. — Optical CM diagram showing Ha emission line stars 
(orange crosses), Class II- type objects (blue circles), and H — K > 
0.6 so urce (red circle). T he thick solid line is the 2 Myr isochrone 
by Gir ardi et al.| ( |2002^ , whi le the thinner so lid lines are PMS 
isochrones tor ditterent ages by |Siess et al.| | 2000||. Das hed lines are 
evolutionary tracks for different masses •^ ISiess et ari (2000>. The 
isochrones have been corrected for distance (i.i kpcj and reddening 
(Ay = 2.70 mag). The arrow shows the reddening vector for Ay 
= 3 mag. 

evol utionary tracks for masses ranging from 0.3 to 3 Mq 
from [Siess et al.| ( |2000[ ) . 

We can see that most of the YSOs lie in the age range 
of 0.5 to 2 Myr, with the mean age ~ 1 Myr. The mass 
range mostly varies from 0.3 to ~ 2 Mq, with a few 
outliers mainly being the Ha emission line sources. The 
PMS star to the right of 0.1 Myr isochrone is likely to 
be a highly extincted source. It should be noted that 
variable extinction along the line of sight and variability 
could be the reason behind the age spread in our figure. 
Similar age spreads have been noted in other s tar- forming 
regions ( |Jose et al.||2008| [Sharma et al.|[2007[ ). 



2.70 mag; see Section 5.3). The figure also shows the 



5.5. Mass spectrum of the stellar sources 

The excess if -band emission can make the YSOs ap- 
pear luminous and hence the K/H — K CM diagram 
is not appropriate for mass estimation. We therefore 
choose the J/J—H CM diagram to minimize the ef- 
fect of this excess emission which would lead to an over 
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Fig. 10.— J/J-H CM diagram for the YSOs. Blue circles denote 
the Class Il-type sources, while red circles denote those with H — 
K > 0.6. The Ha emission line stars have been represented by 
orange crosses. The central ionizing star has been marked by a 
green square . The PMS isochrone for 1 Myr (solid line) from |Siess| 
let al.l (!2000| has been shown. The isochrone was corrected tor a 
distance ot 1.1 kpc. The slanted lines show the reddening vectors 
for different masses. 



estimation of mass. Figure [TO] shows the required CM 
diagram, with the Class Il-like sources detected from the 
CC diagram (in blue circles), sources with H — K > 0.6 
(in red circles), and sources with Ha emission line (in 
orange crosses). The ionizing star has also been shown 
with a green square. The larger statistics here as opposed 
to the optical CM diagram is due to the young nature 
of the sources, which emit at longer wavelengths. From 
the optical CM diagram, we found that the average age 
of the YSOs was about 1 Myr. Also, the dynamical age 
of this H II regio n from radio results is 1.07 ± 0.31 Myr 



(see Section |4.2[ ). Hen ce, we use the iso chrone for 1 Myr 
(black solid line) from Siess et al. (2000) for our analysis. 
The slanting dotted lines show the reddening vectors for 
different masses. 

We can see from the 1 Myr reddening vectors, that 
most of the Ha emission line objects and Class Il-like 
sources lie between 0.1 to 2 Mq, though a few outliers 
seem to be much more massive, of the order of ~ 3 M0. 
Most H — K > sources have masses less than 2 M0. 
It is noteworthy that our mass range is consistent with 
the mas s ra nge derived from the optical CM diagram 
(Section |5.4| , with the Ha emission line stars being to- 
wards the massive end. As can clearly be seen, the YSOs 
with the assumed age, have widely varying colors. This 
is probably an indication of variable extinction, weak 



contribution of excess emission in J and H bands, and 
sources being in different evolutionary stages. It should 
be noted that mass estimations from infrared CM dia- 



grams are prone to systematic errors (Hillenbrand et al. 
20081) due to possible use of different PMS evolution- 



ary models, and reliance on uncertain age and distance. 
Another error source could be the presence of binaries, 
which will brighten a star - thus changing the mass and 
age estimates. 

5.6. Spectral energy distribution 

So far, we have identified the YSOs with the help of 
NIR CC and CM diagrams, and the Ha emission line 
stars. To get an idea about the physical parameters of 
these sources, we modelle d the SEDs using the grid o f 
models and fitting tools of [Robitaille e t al.| (|2006| |2007[). 



The basic model assumes three components of any 
source : a PMS central star, a surrounding flared ac- 
cretion disk, and a rotationally flattened envelope with 
cavities carved out by a bipolar outflow. Different models 
are computed for a suitably large parameter space - con- 
taining manifold combinations of the above mentioned 
compone nts - using a M onte Carlo based radiation trans- 
fer code (I Whitney et al.! 2003a b). SED fitting for differ- 



ent sources is not unique and the number of models for 
each source can only be constrained by having more data 
points from multiwavelength observations. To reduce er- 
ror due to the uncertainties which result from very few 
data points, we modelled only those sources for which 
we had WISE data and at least five data points. The in- 
put parameters were distance range (taken to be 0.9—1.3 
pc), and Ay (taken as 1.1—25 mag). The low er limit was 
chosen to b e 1.1 mag (from the calculation of Bonatto & 
Bica| ( |2009| ) for BDSB 96) as the YSOs should have an ex- 
tinction which is more than the visual absorption for the 
region. The upper limit was chosen as 25 mag because 
the highest individual Ay values from col or-excess cal- 
culations using NICE method (Se ction 1 5. 31) were around 
this value. The SED fitting tool of lRobitallle et al.| ( |2007 ) 
gives a set of well constrained values tor each parameter, 
rather than the values for o nly the best fit for SED. Fol- 
lowing a similar method as Robitaille et al. (2007), we 
considered only those fits to constrain the physical pa- 
rameters for which we had 

- Xrmn < 3 (per data point), (6) 

where is m easure of goodness of fit for each model 
(see Figure [TTj). The final values were subsequently ob- 
tained using the models which satisfied Equation [6] and 
finding a weighted mean value for each parameter, with 
the weight being the inverse square of for each model 
respectively. For some of the sources, the errors for a few 
parameters are large because we are dealing with a large 
parameter space, while we have only a few data points to 
fit the models. The results from our SED modelling are 
given in Table [2j The physical parameters listed are stel- 
lar mass {Mq), stellar age (i*), temperature (T), disk ac- 
cretion rate (Mdisk), foreground visual extinction (Ay), 
and the Xmin P^^ data, point of the be st fit. 

Two of the PMS sources detected by |Linz et aL| ( |2010| ) 
- LDN1657A-2, & LDN1657A-3 - had WISE counter- 
parts. LDN1657A-3 was also detected at NIR wave- 
lengths. Using the flux values given for these sources at 
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TABLE 3 

SED RESULTS FOR, THE YSOS 



ID^ 


RA 
(J2000) 


Dec. 
(J2000) 


Mass log t, log T 
(Mo) (yr) (K) 


log Mdisfe 
(Mo yr-i) 


Av 
(mag) 


^2 ■ 


Sources from Linz et al. (2010) 


1 

2 


07:04:58.030 
07:05:00.665 


-12:16:50.80 
-12:16:45.13 


0.21 ± 0.39 3.37 ± 0.46 3.46 ± 0.07 
1.27 ± 1.49 5.32 ± 0.39 3.61 ± 0.10 


-6.82 ± 0.55 
-7.71 ± 1.85 


16.23 ± 7.82 
10.56 ± 4.20 


4.36 
1.89 


Sources with H — K > 0.6 


19 
45 
49 


07:05:00.919 
07:05:10.395 
07:05:13.678 


-12:15:42.37 
-12:19:25.70 
-12:19:51.96 


2.44 ± 1.00 5.43 ± 0.44 3.67 ± 0.05 
2.91 ± 1.38 4.55 ± 1.27 3.81 ± 0.27 
1.14 ± 0.72 5.05 ± 0.38 3.61 ± 0.03 


-9.15 ± 1.97 
-6.35 ± 1.45 
-8.18 ± 1.11 


7.27 ± 1.14 

1.89 ± 0.78 

2.90 ± 1.83 


1.40 
8.28 
1.20 


Class Il-type sources 



53 


07:04:56.673 


-12:16:12.77 


0.97 ± 0.82 


5.21 ± 0.37 


3.60 ± 0.03 


-8.74 ± 1.20 


2.22 ± 1.16 


3.23 


54 


07:04:57.964 


-12 


17:13.85 


0.19 


± 


0.22 


5.90 


± 


0.48 


3.49 


± 


0.03 


-10.05 ± 1.46 


3.94 ± 0.88 


1.45 


58 


07:05:01.230 


-12 


19:03.69 


1.28 


± 


0.90 


5.15 


± 


0.51 


3.62 


± 


0.04 


-9.01 ± 1.41 


2.98 ± 1.32 


3.34 


64 


07:05:08.269 


-12 


19:15.83 


3.62 


± 


2.24 


5.37 




1.20 


3.87 




0.30 


-8.82 ± 2.73 


15.93 ± 5.87 


0.34 


66 


07:05:09.197 


-12 


19:00.00 


0.86 


± 


0.74 


4.63 


± 


0.72 


3.59 


± 


0.09 


-7.88 ± 1.20 


4.72 ± 2.26 


1.64 


69 


07:05:11.424 


-12 


19:24.36 


0.86 




0.71 


4.72 




0.65 


3.59 




0.08 


-7.98 ± 1.15 


4.45 ± 1.95 


1.50 


70 


07:05:11.541 


-12 


18:11.39 


1.32 




1.35 


4.43 




0.79 


3.60 




0.06 


-7.93 ± 1.12 


3.07 ± 1.37 


2.21 


72 


07:05:11.774 


-12 


16:24.01 


1.00 


± 


0.80 


5.24 


± 


0.58 


3.64 




0.13 


-8.56 ± 1.45 


2.71 ± 1.52 


2.13 


73 


07:05:14.495 


-12 


20:47.57 


0.53 




0.79 


4.93 




0.89 


3.52 




0.14 


-8.76 ± 1.21 


2.45 ± 1.19 


0.03 


75 


07:05:15.324 


-12 


20:02.05 


0.84 




1.09 


4.42 




0.90 


3.58 




0.15 


-7.45 ± 1.23 


2.60 ± 1.92 


2.15 


76 


07:05:17.419 


-12 


18:58.69 


0.81 


± 


1.19 


4.85 


± 


1.15 


3.56 


± 


0.22 


-8.20 ± 1.42 


2.21 ± 1.01 


0.04 



Ha emission line sources 



"79" 
80 
82 
86 
87 
88 
89 
90 
91 
92 
93 



07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 
07:05: 



06.330 
09.652 
12.292 
17.904 
18.494 
19.389 
19.592 
20.068 
20.728 
21.912 
22.359 



-12:15: 
-12:19: 
-12:18: 
-12:15: 
-12:15: 
-12:20: 
-12:17: 
-12:19: 
-12:21: 
-12:19: 
-12:17: 



36.38 
56.32 
38.27 
40.26 
09.22 
28.71 
31.48 
12.59 
42.69 
07.19 
18.63 



2.27 ± 0.67 
0.89 ± 1.18 
0.94 ± 0.92 
4.60 ± 1.03 
2.30 ± 0.72 
2.51 ± 0.81 
2.33 ± 0.86 
3.93 ± 0.15 
2.64 ± 0.78 
0.26 ± 0.18 
1.07 ± 0.88 



6.63 ± 0.36 

4.60 ± 0.82 
5.68 ± 0.30 

5.61 ± 0.07 
6.04 ± 0.58 
6.55 ± 0.31 
5.58 ± 0.42 
6.63 ± 0.13 
5.57 ± 0.31 
5.83 ± 0.29 
6.14 ± 0.45 



3.85 ± 0.15 
3.58 ± 0.06 

3.61 ± 0.04 

3.78 ± 0.09 

3.79 ± 0.17 
3.82 ± 0.16 

3.67 ± 0.05 
4.15 ± 0.01 

3.68 ± 0.05 
3.51 ± 0.02 

3.62 ± 0.06 



-9.95 ± 1.95 
-7.90 ± 0.98 
-8.76 ± 1.15 
-7.48 ± 1.00 
-9.08 ± 1.49 
-9.41 ± 2.01 
-9.61 ± 1.32 
-12.90 ± 0.60 
-9.15 ± 1.64 
-10.70 ± 1.64 
-9.20 ± 1.18 



5.11 ± 1.05 
3.24 ± 1.59 
1.46 ± 0.52 
2.01 ± 0.60 
2.31 ± 1.05 
2.23 ± 0.92 

2.12 ± 0.58 
3.83 ± 0.21 
1.91 ± 0.66 
1.34 ± 0.34 
2.06 ± 1.00 



1.15 
1.52 
2.62 
0.54 
0.24 
2.68 
1.27 
1.04 
2.61 
1.25 



Ho emission line sources outside JNIK i^bV 


96 


07:05:23.669 -12:22:36.08 


1.69 ± 0.64 6.22 ± 0.38 3.71 ± 0.13 


-9.96 ± 1.58 


2.46 ± 1.11 


1.64 


^ The ID 


no. refers to the ID in Table[2] 
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Fig. 11.— SED fitting using the model of |Robitaille et al.| ( [2007| l for (left) LDN1657A-3, from |Linz et al.| ||2010'), and [right) a typical 
Class Il-type source. The ID no. corresponds to tne IJJs m tne YbU table (Table[2]l. The black dots are our tne oata points. The solid black 
line shows the best fitted model, while the grey lines show the consecutive good fits for — Xmin {P^^ data point) < 3. The dashed line 
shows the fit for the photosphere of the central source for the best fit model. The photosphere fit is assuming that there is no circumstellar 
dust, however, interstellar reddening has been accounted for. 



longer wavelengths from Linz et al. (20101, we could fur- 
ther constrain the physical parameters or these sources. 
The source with ID "1" in our SED table (Table [2| 
gives the values for LDN1657A-2, while "2" gives those 
for LDN1657A-3. Figure [ll] shows examples of SEDs 
with the resulting models for LDN1657A-3 and a typical 
Class-II like source. 
It must be kept in mind that the SED results are only 



representative of the actual values and should not be 
taken too literally. Out of a database of 200000 YSO 
models (each for a different set of parameter values), 
the model simply finds the closest matches for a par- 
ticular source entered by the user. Any grid of models 
will have inherent limitations, which could give rise to 
pseudo-trends in the results. The models are for isolated 
objects and thus could be misleading when we have un- 
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resolved multiple point sources. Moreover, the model as- 
sumes same phys ics for stellar masses from 0.1 to 50 Mq. 
Robitaille (2008) provides a succinct summary of the 



caveats inherent in this modelling. Nevertheless, we can 
still use the results to glean general information about 
the mass and age ranges of our sources. As can be seen 
from Table [2] most of the sources have stellar ages in the 
range 0.05 to 1 Myr and masses in 0.2 to 2.5 Mq range. 
The outliers are expectedly the more evolved Ha emis- 
sion line stars, whose ages vary upto 5 Myr, and stellar 
masses are ~ 2.5—5 Mq. We see a few very low mass 
(3 sources with masses < 0.5 Mq) and extremely young 
sources (9 sources with ages < 10^ yr, out of which 6 are 
CTTS) as well. The disk accretion rate for most of the 
Ha emission line sources is in the range of ^ 10~^-10~^^ 
M0 yr-i, which is less than that of CTTS - IQ-'^-lQ-^ 
Mq yr"-'^, though it would be prudent to treat these val- 
ues with caution owing to a lack of FIR and submillime- 
tre data points. A few of our sources, ID no.s 39, 55, 
61, 81, 83, 94, 95, and 97 (Table [2]), were better fitted 
by stellar models of a heavily extincted central star, and 
are thus not included here. Hence, we find that our age 
and mass estimations are mostly consistent with those 
derived from the optical CM diagram and the mass spec- 
trum above. 



5.7. Spatial distribution of the YSOs 

Figure [12] shows the spatial distribution of the YSOs 
from our IRSF and grism slitless spectroscopy data. Red 
circles denote the sources with H — K > 1, sources with 
1 > H — K > 0.6 have been shown with magenta crosses, 
blue circles denote the Class Il-like objects detected from 
our CC diagram, while the white crosses show the loca- 
tion of the Ha emission line stars from grism slitless spec- 
troscopy. The central ionizing source HP 536 23 is shown 



with a green square. Forbrich et al. (2004) had found 



a young, intermediate mass protostar named UYSOl at 
the north-west border of the H ii region, towards the cold 
dark cloud LDN1657A. The location of UYSOl, which 



was f urther resolved into two protostars (Forbrich et al. 
2009|, has been denoted by a green plus sign. In addi- 



tion, a recent Herschel PACS and SPIRE mapping of the 
region around UYSOl identified five cool and compact 
FIR sources, of masses of t he order of a few solar masses, 
hidden in the dark cloud fLinz et al. 2010). They have 
been marked by green circles. Most of the YSOs iden- 
tified by us lie in the molecular cloud reg ion traced by 
^^CO molecular emission (see Section 4.3 & Figure [6|i, 
along the concave arc facing the north-eastern corner. 

We have divided the region into three sub-regions 
which show increasing infrared excess from the south- 
east towards the north-west direction. This indicates 
an evolutionary sequence, thereby suggesting a possi- 
ble decreasing age gradient towards that direction. The 
youngest stars, i.e. those with H — K> 0.6 mostly lie to- 
wards the north-west direction, followed by a concentra- 
tion of Class Il-like objects and then the Ha emission line 
stars towards the ionizing region. Even among the stars 
with H — K > 0.6, the redder sources {H — K > 1) were 
found to be lying towards the north-west, with a decreas- 
ing color sequence towards the ionizing source. Also, the 
5 sources found by Hersc hel mapping are cooler than the 
massive object UYSOl (Linz et al. 2010), thereby indi- 




RA (J2000) 



Fig. 12. — Color composite showing the spatial distribution of 
the YSOs, made using IRSF J (blue), H (green), and Kg (red) 
band images. The Ha emission line stars have been denoted by 
white crosses, Class Il-like sources by blue circles, sources with 
1 > H — K > 0.6 with magenta crosses, and sources with H ~ 
K > 1 hy red circles. The green square marks the position of 
the central ionizing star HD 53623. UYSOl (Forbrich et al. 2004| 
has been ma rked with a green plus sign, while the cool, compact 
sources from |Linz et al.| | |2010l l by green circles. 1280 MHz GMRT 

contours (~ 25 resolution) showing the champagne flow have been 
shown in dashed white lines. The three sub-regions denoted by the 
yellow boxes show an increasing fraction of redder sources in the 
sub-regions as we move from the ionizing star towards the north- 
western direction. The dashed yellow lines are drawn at a distance 
of 80 and 220 from the ionizing star, respectively. 

eating that they are younger than UYSOl. We tried to 
quantify the infrared excess {H — K) for increasing dis- 
tances of the sources from the ionizing star HD 53623. 
Figure 13 shows a plot oi H — K color versus distance. 
The graph is divided into three regions as shown by three 
boxes in Figure 12, with the average H — K color calcu- 
lated for each of them. The average color was found to 
increase as we move away from the ionizing star, thereby 
indicating that the sources in general, are younger at 
larger distances from the central region. 

This distribution of YSOs, with relatively larger NIR 
excess at larger radial distances, could be indicative of 
probable triggered star formation, which seems to have 
propagated from the ionizing source towards the north- 
west. 

6. STAR FORMATION SCENARIO 

Star formation in a molecular cloud is more active dur- 
ing the first few million years of cloud time. Subsequent 
to that, the gas in molecular cloud dissipates and further 
star formation no longer takes place. However, these 
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Fig. 13. — The variation oi H — K color with distance from the 
star HD 53623. The Ha emission line stars have been denoted 
by black circles, Class Il-like sources by blue circles, sources with 
1 > H—K > 0.6 with magenta crosses, and sources with H—K > 1 
by red circles. Green circle denotes LDN1657A-3 source. 

types of molecular cloud complexes can still have star 
formation at the interface of the H ii region created by 
OB stars and the molecular gas. In Sh2-297 region, the 
spatial distributio n of the YSOs and infrared excess cal- 
culations (Section |5.7| suggest an evolutionary sequence 
of sources from the ionizing region towards the cold, dark 
cloud LDN1657A, with a majority of YSOs distributed 
at the western and north-western border of Sh2-297. 

One of the conditions for the massive star HD 53623 
to have influenced the formation of YSOs in the region 
is that the age of the H ii region and associated sources 
should be of the order of, or greater than, the age of the 
YSOs in the triggered region. However, it is difficult to 
estimate the age of YSOs as very few of them have op- 
tical counterparts. Nonetheless, we tried to estimate the 
age of the few sources which have optical counterparts 
in Section 5.4 Age estimation using IR data could be 
misleading because IR emission arises not only from the 
photosphere, but also the circumstellar disk of the YSOs. 
We tried to constrain the age of embedded YSOs with the 
SED fitting models, though due to the lack of sufficient 
data points, the estimated ages are uncertain. However, 
they can still be used as quantitative indicators of stellar 
youth. Our analysis shows the mean age of YSOs to be ^ 
1 Myr, which is of the order of the dynamical age of the 
H II region of 1.07 ± 0.31 Myr, considering the uncertain- 
ties in age estimation. Our SED results (Table [2| also 
point out that barring 6 sources (all being Ha emission 
line sources; possibly associated with the evolved H ii re- 



gion), the other 22 sources have ages less than 1.07 Myr 
(majority of them are located in the dark cloud at the 
western b order of Sh2-297). Th e protostellar source de- 



tected by [Forbric h et al.| (|2009 ) 
sources detected 



' Lmz et al 
star formation occurring m the 



, and the cool, compact 
|2010) also show recent 
'dark cloud. The dynami- 



cal age of the outf l ow for UYSOl, which was constrained 



by Forbrich et al 
the immediate in 



( |2009[ ) as < 10* yr, and its location at 
erface of Sh2-297 supports our notion 
of triggered star formation by the expanding H il region. 

To examine if the collect-and-coUapse due to the ioniz- 
ing star is possibly at work in the region, we compare 
the observed propertie s with the analytical model by 



Whitworth et al. ( 1994 Section 5) for star formation at 
the periphery ot an expanding H ii region. The model 
gives the expressions for the time when the fragmentation 
starts (t frag) 7 Slid the radius at the time of fragmenta- 
tion (ifrag)- The valuc of Os (isothermal sound speed in 
the shocked layer), which varies from 0.2 to 0.6 km s^^, 
was taken to be 0.2 km s^^ for our calculations. Using 
the valu es o f and Uq calculated in previous section 
(Section 4.2 1, we obtained i/rag ~ 1-27 Myr and rfrag ~ 
1.70 pc. I'liis would suggest that fragmentation would 
start at a radius of 1.70 pc, whereas the stellar sources 
found are within a 1.70 pc radius of the massive star HD 
53623. However, the distance of the YSOs is the pro- 
jected distance and the actual distance could be larger. 
Taking the present mean age of the YSOs as ~ 1 Myr 
and the fragmentation time as calculated, it would imply 
that the age difference between the YSOs and HD 53623 
(-- 5 ± 3 Myr) is > 2 Myr. Though this tfrag calcula- 
tion suggests that this H il region could be old enough 
to lead to star formation by collect-and-collapse process, 
the morphology of the region makes this process unlikely. 
Moreover, we have assumed the lowest value of isother- 
mal sound speed for our calculations (0.2 km s~^). A 
higher value of Os will only further increase the fragmen- 
tation time and radius and make a collect-and-collapse 
scenario unlikelier. Also, the YSOs are not exclusively 
associated with the PAH emission that we see around 
this region (Figure [6| . PAH regions represent the com- 
pressed neutral matter between the ionization and shock 
fronts. The YSOs extend well beyond the PAH emission, 
as is evident from our MSX 8.28 /im contours (which, in 
addition, contains continuum emission too). This is dis- 
tinctly separate from the star forming regions of, for ex- 
ample, Deharveng et al. (2005) where young sources and 
condensations have been observed along the PAH ring. 
It should be noted that, in an ideal case, one would ex- 
pect regular ly-spaced-clumps/protostars of very similar 
ages distributed along the periphery of the H ii region, 
mostly along the material swept-up by the expanding 
ionization front. This is definitely not the case here. It 
therefore appears that the collect-and-collapse process 
due to the star HD 53623 is most likely not involved in 
the star formation in this region. 

Alternatively, it is possible that the star formation ob- 
served at the north-western border of Sh2-297 is proba- 
bly due to some other processes. The formation of YSOs 
in the dark cloud area could partly be due to the influ- 
ence of expanding H ii region Sh2-297 by RDI (radiation 
driven implosion) mechanism (most likely upto the edg e 
and its vicinity) (Bertoldi| 1989} [Lefloch fc Lazareff 1 1 994 j ) , 
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and partly due to supernova induced star formation as 
discussed by Herbst & Assousa ( 1977) (the distant YSOs 
possibly having tormed by this process). The YSOs are 



distributed along the concave arc (see Section 5.7) of the 
molecular cloud traced by ^^CO, and seem to have an 
age sequence. It appears that the supernova might have 
swept the molecular material and led to its piling up 
along the arc. Therefore, as the material might have col- 
lected along the arc, beginning from the edge of Sh2-297 
towards the north-west, star formation could have oc- 
cured at different epochs. The plausible lo cation of the 



supernova (a igso ~ 07 08™, 



950 



-11.2°, 



Reynolds & 



Ogden [l978 1 a nd the derived age of the supernova rem- 
03TMyr ( [Herbst fc Assousa|1977| - lend weight to 



nant 

this hypothesis.' 

To better understand the star formation and the ef- 
fect of the massive star on the Sh2-297 region, we need 
much deeper, wider and longer wavelength observations 
for better statistics (so as to constrain mass and lumi- 
nosity functions of the embedded YSOs), and to carry 
out a study of YSOs along the arc of molecular material 
connecting Sh2-297 to the larger region. It would also 
be helpful to carry out a radial velocity study and spec- 
troscopic observations for a wider area to constrain the 
membership, age and spectral type of the YSOs. This 
would confirm their relative evolutionary status', iden- 
tify possible clusters, and help put the relevant triggered 
star formation process on a firm footing. 

7. CONCLUSIONS 

In this paper, we have carried out a multiwavelength 
study of the star-forming region Sh2-297. Our main con- 
clusions are as follows : 

1. We identified the spectral class of the central ion- 
izing star HD 53623, using optical spectroscopic 
data. The prominence of He I lines in our spec- 
tra indicate a spectral type of BOV. Our optical 
spectroscopic results are consistent with the radio 
continuum data results. 

2. Using the radio data for 610, 1280, & 1420 MHz, 
the physical parameters of the region were derived. 
Taking the extent of the H ii region to be 1.6 pc 
(5 ), the emission measure was calculated to be 
9.15 ± 0.56 X 10^ cm~^ pc, and thereby the elec- 
tron density as 756 ± 46 cm~'^. These low val- 
ues of emission measure and the electron density 
confirm that this is indeed a classical H ii region. 
The stromgren radius was calculated to be 0.051 ± 
0.002 pc and the dynamical timescale to be 1.07 ± 
0.31 Myr. 

3. The grism slitless spectroscopic image was used to 
identify the Ha emission line stars; NIR CC dia- 



gram [J-H/H-K) to identify the Class II- type ob- 
jects, and the NIR CM diagram (K/H-K) to iden- 
tify the red sources {H — K > 0.6) in our field. The 
average extinction in the field was calculated to be 
2.70 mag. 

4. We estimated the age and mass of the YSOs using 
our optical CM diagram ( V/V-I), besides using the 
NIR mass spectrum to estimate the mass range. 
Age was found to be in the range 0.5 to 2 Myr, 
with the average age being ^ 1 Myr. The mass 
was estimated to be ~ 0.1 to 2 Mq, with a few 
outliers ^ 3 Mq. 

5. We studied t he evolutionary st atus of some of the 
sources from Linz et al. (2010), red sources. Class 
Il-type objects, an d Ma emission line star s using 
the SED model of [Robitaille et al.| ( |2007[ ). The 
models predict - for most of the sources - stellar age 
in the range 0.05—1 Myr, stellar mass in the range 
0.2-2.5 Mq, and disk accretion rate ~ 10"'^-10"" 
MQyr~^. A few Ha emission line stars, which are 
most evolved YSOs in the field, have stellar ages 
upto 5 Myr, and masses ~ 2.5—5 M0. 

6. The spatial distribution of the YSOs indicates a 
possible evolutionary sequence, in the sense that 
youngest sources are distributed away from the ion- 
izing source. This could be evidence in support of 
triggered star formation in the region. The star 
formation seems to have propagated from the ion- 
izing source towards the cloud LDN1657A, which is 
further supported by the prese nce of massive pro - 
tostar UYSOl discovered by Forbrich et a ll (!2004L 
and co ol, compact sources discovered by 'Linz ct m] 
(2010), which are located away from the ionizing 



source. 
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Mallick et al. 



TABLE 2 

Optical &; infrared magnitudes for the YSOs 



ID 



RA 
(J2000) 



Dec. 
(J2000) 



H 



wl 



w2 



LDN1657A-2 fc 3 from Linz et al. (2010) 



07:04:58.030 
07:05:00.665 



-12:16:50.80 
-12:16:45.13 



14.12 =t 0.08 10.78 ± 0.02 S: 
11.12 ± 0.01 9.73 ± 0.01 7. 



16.85 ± 0.10 14.67 ± 0.05 12.96 ± 0.04 



YSOs with - K > 0.6 ~ 
16.62 ± 

16.82 ± 
17.32 ± 
16.67 
14.95 
15.87 
13.32 
17.14 
16.60 
16.27 ± 
16.73 ± 
16.12 ± 
15.69 ± 
16.25 ± 
16.80 ± 
16.86 ± 

10.83 ± 
16.83 
16.95 
16.99 
16.29 
17.34 
16.09 
15.10 ± 
15.98 ± 
16.88 ± 
17.15 ± 
16.91 ± 
16.55 ± 
16.46 ± 
16.85 ± 
15.48 ± 



3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 



07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:04; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 



55.818 
56.305 

56.334 
56.580 
56.634 
56.898 
57.702 
57.927 
58.407 
58.749 
59.150 
59.994 
00.319 
00.599 
00.826 
00.881 
00.919 
01.313 
01.873 
02.311 
02.637 
02.833 
03.666 
03.832 
04.448 
05.984 
05.995 
06.044 
06.075 
06.216 
06.735 
07.139 
08.827 
08.890 
09.181 
09.283 
09.326 
09.749 
09.811 
10.071 
10.124 
10.358 
10.395 
10.481 
11.475 
11.946 
13.678 
17.596 
17.668 
23.035 



-12:17; 
-12:17; 

-12:16; 
-12:15; 
-12:16; 
-12:16; 
-12:16; 
-12:20; 
-12:17; 
-12:16; 
-12:15; 
-12:15; 
-12:16; 
-12:18; 
-12:18; 
-12:18; 
-12:15; 
-12:17; 
-12:17; 
-12:17; 
-12:16; 
-12:16; 
-12:18; 
-12:18; 
-12:17; 
-12:17; 
-12:16; 
-12:17; 
-12:18; 
-12:17; 
-12:17; 
-12:19; 
-12:17; 
-12:15; 
-12:19; 
-12:17; 
-12:18; 
-12:20; 
-12:19; 
-12:18; 
-12:18; 
-12:18; 
-12:19; 
-12:18; 
-12:19; 
-12:18; 
-12:19; 
-12:20; 
-12:20; 
-12:19; 



37.90 
24.14 
14.28 
14.47 
59.69 
39.00 
20.06 
36.37 
47.23 
18.35 
17.78 
41.34 
43.75 
52.50 
46.47 
26.03 
42.37 
22.48 
11.83 
17.60 
44.26 
12.72 
24.20 
00.40 
09.32 
21.69 
45.27 
17.39 
35.45 
33.02 
52.68 
23.87 
54.66 
34.44 
25.46 
24.01 
42.24 
02.05 
26.43 
51.91 
29.40 
09.88 
25.70 
30.70 
16.58 
48.14 
51.96 
23.48 
26.39 
56.15 



17.39 



17.65 
17.56 



± 0.04 



± 0.03 
± 0.03 



19.65 



0.01 15.29 



± 



17.55 ± 0.05 
0.01 12.47 ± 0.01 



16.64 



± 0.02 



18.25 ± 0.07 



17.59 ± 0.05 

17.69 ± 0.03 

16.00 ± 0.01 

17.70 ± 0.04 
14.38 ± 0.04 
17.45 ± 0.04 



16.65 ± 0.01 
17.46 ± 0.03 



15.46 ± 0.01 
15.89 ± 0.01 

17.47 ± 0.04 
17.03 ± 0.02 



16.31 ± 



16.i 



± 



14.15 ± 
16.31 ± 
12.84 ± 
16.22 ± 
16.97 ± 
18.02 ± 
14.53 ± 
16.26 ± 
16.64 ± 
17.59 ± 
16.55 ± 
13.99 ± 
14.52 ± 
16.88 ± 
16.66 ± 
16.28 ± 



003" 

0.03 

0.04 

0.02 

0.01 

0.07 

0.01 

0.05 

0.03 

0.03 

0.03 

0.02 

0.03 

0.02 

0.05 

0.04 

0.01 

0.03 

0.04 

0.03 

0.02 

0.04 

0.02 

0.01 

0.01 

0.03 

0.04 

0.04 

0.03 

0.03 

0.03 

0.01 

0.02 

0.06 

0.01 

0.02 

0.06 

0.02 

0.04 

0.10 

0.01 

0.02 

0.06 

0.07 

0.06 

0.01 

0.01 

0.04 

0.03 

0.02 



15.51 
15.47 

16.18 
15.80 
13.60 
15.09 
12.21 
16.12 
15.27 
15.15 
16.05 
15.34 
13.44 
15.52 
16.12 
16.10 
10.07 
15.45 
15.87 
16.21 
14.99 
15.72 
14.27 
14.35 
14.31 
15.37 
15.39 
15.65 
14.59 
15.48 
15.85 
14.34 
15.60 
16.00 
13.26 
15.63 
12.05 
15.52 
16.06 
16.19 
13.40 
15.61 
15.66 
15, 
15.38 
13.32 
13.83 
16.14 
15.86 
15.64 



± 0.09 

■ ± 0.07 

± 0.13 
0.10 
0.02 
0.06 
0.01 
0.13 
0.06 
± 0.06 
± 0.14 

: ± 0.06 

0.02 
± 0.08 
± 0.13 

I ± 0.14 
0.01 
0.08 
0.12 
0.13 
0.05 
0.09 
0.03 
± 0.03 
± 0.03 

■ ± 0.06 
i ± 0.07 
I ± 0.09 
i ± 0.03 
; ± 0.08 

± 0.11 
0.03 
0.08 
0.13 
0.01 
± 0.09 
± 0.05 
± 0.08 
i ± 0.14 
i ± 0.14 
I ± 0.02 
± 0.09 
1 ± 0.10 
± 0.10 
0.13 
0.01 
0.02 
0.13 
0.11 



14.78 ± 0.11 13.69 ± 0.10 

14.63 ± 0.08 

12.64 ± 0.03 

13.71 ± 0.06 12.49 ± 0.08 

11.06 ± 0.02 10.91 ± 0.02 

14.78 ± 0.10 14.21 ± 0.16 



14.59 ± 0.10 
9.61 ± 0.01 
14.72 ± 0.10 
14.82 ± 0.16 



9.47 



± 



0.01 



14.46 ± 0.18 
13.23 ± 0.05 
13.12 ± 0.06 

13.26 ± 0.05 12.57 ± 0.05 



14.33 ± 0.24 
13.69 ± 0.10 
12.15 ± 0.05 



14.65 
13.27 

15.72 
10.1 

11.20 
11.81 



± 0.08 



± 0.14 

± 0.09 12.06 

± 0.27 

± 0.06 10.12 

± 0.10 10.85 

± 0.13 11.10 



± 0.07 



± 0.06 



± 0.12 5. 



± 0.10 



Class 11-typc sources 



"53" 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 



07:04; 
07:04; 
07:04; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 



56.673 
57.964 
59.938 
00.844 
01.117 
01.230 
02.673 
03.364 
04.206 
04.644 
05.166 
08.269 
08.771 
09.197 
10.638 
10.818 
11.424 
11.541 
11.642 
11.774 



-12:16; 
-12:17: 
-12:17; 
-12:20; 
-12:16; 
-12:19; 
-12:15; 
-12:15; 
-12:18; 
-12:19; 
-12:18; 
-12:19; 
-12:18; 
-12:19; 
-12:19: 
-12:19; 
-12:19; 
-12:18; 
-12:19: 
-12:16; 



12.77 
13.85 
24.19 
06.48 
50.42 
03.69 
44.75 
53.58 
45.13 
10.39 
19.49 
15.83 
51.26 
00.00 
15.02 
43.90 
24.36 
11.39 
30.95 
24.01 



20.19 ± 0.03 16.94 ± 0.01 



15.81 
16.08 
17.61 
17.20 
17.35 
14.83 
16.70 
17.78 
17.90 
15.96 
17.43 
16.16 
16.67 
17.09 
17.04 
17.16 
16.93 
14.63 
15.10 
15.34 



JHTT 
± 0.01 
± 0.05 
± 0.04 
± 0.04 
± 0.01 
0.02 
0.04 
0.06 
0.02 
0.04 
0.01 
0.02 
± 0.03 
± 0.02 
± 0.04 
± 0.04 
± 0.01 
± 0.01 
± 0.01 



14.25 ± 

15.07 ± 
15.61 ± 
16.32 ± 

16.08 ± 
13.60 ± 
15.84 
16.75 
15.30 
14.82 
15.88 
13.69 
15.41 ± 
15.29 ± 
15.89 ± 
15.83 ± 
15.13 ± 

13.09 ± 
13.96 ± 
13.72 ± 



0.01 
0.03 
0.02 
0.03 
0.01 
0.01 
0.03 
0.01 
0.01 
0.02 
0.01 
0.02 
0.01 
0.03 
0.02 
0.01 
0.01 
0.01 
0.01 



13.29 
14.37 
14.43 
15.70 
15.13 
12.85 
15.28 
16.09 
13.80 
14.13 
14.82 
12.16 
14.64 
13.95 
15.10 
14.89 
13.76 
12.15 
13.14 
12.52 



± 0.02 
± 0.04 
± 0.04 
± 0.09 
± 0.06 
± 0.01 
± 0.07 
± 0.14 
± 0.02 
± 0.02 
± 0.04 



12.04 it 0.02 11.43 ± 0.02 
14.14 ± 0.06 12.92 ± 0.06 
13.17 ± 0.04 12.94 ± 0.07 



20.20 



± 0.03 16.73 



± 0.01 



11.77 ± 0.02 11.11 



20.02 ± 0.05 16.58 ± 0.03 



20.78 ± 0.05 17.95 ± 0.02 



0.01 
0.04 
± 0.02 
I ± 0.06 
i ± 0.05 
1 ± 0.04 
± 0.01 
: ± 0.01 
± 0.01 



12.44 ± 0.04 11.37 ± 0.03 
13.46 ± 0.13 

10.41 ± 0.06 9.03 ± 0.03 

11.16 ± 0.08 10.86 ± 0.17 



11.04 ± 0.10 10.45 ± 0.10 
10.04 ± 0.08 

11.27 ± 0.02 10.29 ± 0.02 7. 



± 0.02 
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ID 


RA 


Dec. 


V I J 


H 


Ks 


wl 


w2 




(j2noo) 


(J2000) 












73 


07;05;14..195 


-12:20:47.57 


15.25 ± 0.01 


14.16 ± 0.01 


13.47 ± 0.02 


12.43 ± 0.12 


11.65 ± 0.12 


74 


07;05;14.610 


-12:19:42.29 


15.07 ± 0.01 


13.62 ± 0.01 


12.74 d= 0.01 






75 


07:05:15.324 


-12:20:02.05 


17.17 ± 0.04 


15.89 ± 0.01 


14.95 ± 0.05 


11.61 ± 0.22 


11.30 ± 0.27 


76 


07:05:17.419 


-12:18:58.69 


14.42 ± 0.01 


13.49 ± 0.01 


12.88 ± 0.01 


11.95 ± 0.04 


11.21 ± 0.06 


77 


07:05:18.113 


-12:19:25.47 


17.25 ± 0.03 


15.90 ± 0.01 


14.94 ± 0.05 






78 


07:05:18.950 


-12:20:13.53 


16.03 ± 0.05 


14.93 ± 0.03 


13.97 ± 0.03 




9.47 ± 0.05 



Ha emission line sources 



"79" 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 



07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 
07:05; 



06.330 
09.652 
10.920 
12.292 
15.359 
15.840 
16.265 
17.904 
18.494 
19.389 
19.592 
20.068 
20.728 
21.912 
22.359 



-12:15; 
-12:19; 
-12:15; 
-12:18; 
-12:19; 
-12:19; 
-12:20; 
-12:15; 
-12:15: 
-12:20: 
-12:17: 
-12:19: 
-12:21: 
-12:19: 
-12:17: 



36.38 
56.32 
22.46 
38.27 
39.02 
45.09 
17.61 
40.26 
09.22 
28.71 
31.48 
12.59 
42.69 
07.19 
18.63 



17.43 ± 0.01 
20.07 ± 0.03 



14.82 
16.18 



± 0.01 
± 0.01 



12.95 ± 0.01 
14.05 ± 0.01 
15.11 ± 0.01 



12.01 ± 0. 
12.77 ± 0. 
14.49 ± 0. 
12.06 ± 0. 
12.27 ± 0. 
12.55 ± 0. 
13.37 ± 0. 
13.65 ± 0. 
12.85 ± 0. 
11.11 ± 0. 
12.51 ± 0. 
10.92 ± 0. 
12.84 ± 0. 
13.44 ± 0. 
13.30 ± 0. 



mr 

01 
01 
01 
06 
01 
01 
02 
01 
01 
01 
01 
01 
01 
01 



12. 
14. 
11. 
11. 
12. 
12. 
13. 
11. 
10. 
12. 
10. 
12. 
13. 
12. 



33 ± 0.01 



10.44 ± 0.01 
10.91 ± 0.07 
13.61 ± 0.04 
11.05 ± 0.11 
10.93 ± 0.12 



12.86 ± 0.03 
10.77 ± 0.01 
10.09 ± 0.06 
11.40 ± 0.01 
9.82 ± 0.01 
11.03 ± 0.02 
12.53 ± 0.03 
12.25 ± 0.02 



9.84 ± 0.01 
10.21 ± 0.07 
13.12 ± 0.05 
10.53 ± 0.12 
10.27 ± 0.10 



12.50 ± 0.04 
9.97 ± 0.01 
9.50 ± 0.05 
10.74 ± 0.02 
9.36 ± 0.01 
10.38 ± 0.02 
12.02 ± 0.03 
11.62 ± 0.03 



0.01 
0.03 



16.74 ± 0.01 14.37 ± 0.01 



15.41 ± 0.14 14.39 ± 0.02 



15.27 
17.35 
14.37 
16.82 
13.86 
17.11 
19.80 
18.23 



± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.02 
± 0.01 



14.17 
14.85 
12.92 
14.59 
12.69 
15.14 
15.92 
15.39 



15.73 
14.14 



± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 
± 0.01 



13.01 
13.13 
13.41 
14.32 
13.89 



0.01 
0.06 
0.01 
0.01 
0.01 



± 

± 
± 
± 

14.00 ± 0.01 
11.83 ± 0.01 
13.38 ± 0.01 
11.67 ± 0.01 
13.83 ± 0.01 
14.30 ± 0.01 
14.17 ± 0.01 



50 ± 0.01 
99 ± 0.05 
10 ± 0.01 



0.01 
0.03 



95 ± 0.01 
79 ± 0.01 
10 ± 0.01 
40 ± 0.01 
15 ± 0.01 
01 ± 0.01 
84 ± 0.01 



emission line sources outside INIH I'oV 

^- 14.15 ± 0.04 13.26 ± 0.04 12.77 ± 0.03 12.26 ± 0.02 11.76 ± 0.02 

± 0.01 14.24 ± 0.03 13.30 ± 0.03 12.77 ± 0.03 12.15 ± 0.02 11.53 ± 0.02 

± 0.01 13.25 ± 0.03 12.48 ± 0.02 12.16 ± 0.02 11.89 ± 0.02 11.46 ± 0.03 

13.98 ± 0.02 13.14 ± 0.02 12.73 ± 0.03 12.04 ± 0.05 11.34 ± 0.07 



95 
96 
97 



07:05:12.139 
07:05:16.154 
07:05:23.669 
07:05:37.457 



-12:14:41.57 
-12:14:55.90 
-12:22:36.08 
-12:21:17.32 



18.57 ± 0.01 
16.18 ± 0.01 



NIR magnitudes taken from 2MASS catalogue 



